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ABSTRACT 

A series of flight tests were conducted on the PC 9/A aircraft, A23-045, 
at the Royal Australian Air Force’s Aircraft Research and Development Unit. 
System identification techniques were applied to the data obtained from these 
flight tests to determine the stability and control derivatives of the aircraft. 
The lateral results for the aircraft in cmise configuration are presented in this 
report and comparisons are made with empirical and ground based estimates. 
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A Correlation between Flight-determined Lateral 
Derivatives and Ground-based Data for the Pilatus PC 9/A 
Training Aircraft in Cruise Configuration 


EXECUTIVE SUMMARY 

Air Operations Division (AOD) has developed, or acquired, a number of fixed-wing 
flight dynamic models and is also responsible for providing advice to the Australian De¬ 
fence Organisation (ADO) on flight simulator flight dynamic model requirements. The 
models generally make use of extensive static and dynamic stability and control derivative 
databases, in addition to engine and flight control models. The static model data may be 
obtained from wind t unn el testing, whilst the dynamic data is traditionally obtained from 
flight tests using system identification techniques. 

The system identification techniques used to estimate aerodynamic derivatives of con¬ 
ventional aircraft are well established. The major requirement of these techniques is high 
fidelity measurements of manoeuvre input (i.e. control surface deflections), and aircraft 
response (i.e. angular rates and linear accelerations), as well as air data including airspeed, 
altitude, angle-of-attack and angle-of-sideslip. 

Following the completion of the AOD PC 9/A wind tunnel tests, the requirement ex¬ 
isted for a dedicated flight dynamic modelling flight test program to both validate the 
flight dynamic model of the PC 9/A and to provide dynamic derivative estimates which 
were unobtainable in the AMRL wind tunnel. An instrumentation suite, including an air 
data probe for the direct and accurate measurement of angle-of-attack, angle-of-sideslip, 
temperature, and static and dynamic pressure, was fitted to the aircraft and a test ma¬ 
noeuvre matrix was designed specifically for the purpose of gathering flight dynamic data. 
This flight test program was conducted at the Aircraft Research and Development Unit 
(ARDU) during 1998/99. 

This report details the analysis of the lateral manoeuvres carried out with the aircraft 
in the cruise configuration. The static and dynamic derivatives thus obtained are com¬ 
pared with wind tunnel estimates as well as a number of empirical estimates obtained 
from alternative sources. A discussion of some of the difficulties encountered during the 
estimation process is also included. The data obtained from these tests will be used in 
the development of the PC 9/A flight dynamic model for the purpose of enhancing AOD 
support for the PC 9/A fleet, including possible upgrades to the part-task tradner. 
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1 Introduction 

The Pilatus PC 9/A is one of a number of high performance turbo-prop aircraft currently 
operated by the Royal Australian Air Force (RAAF). Work currently underway in the Air 
Operations Division (AOD) of the Defence Science and Technology Organisation (DSTO) 
includes investigations into the propeller power effects of such aircraft, and the PC 9/A 
was considerd a suitable platform for study. A six degree-of-freedom flight dynamic model 
of the PC 9/A has been developed by AOD for use in these investigations, as well as for 
pilot-in-the-loop simulations in the Air Operations Simulation Centre and incident and 
accident investigations. 

The development of a flight dynamic model requires information on the static and dynamic 
stability and control derivatives of the aircraft, as well as flight control laws and physical 
properties. The static data for the PC 9/A were collected during both power-off and power- 
on testing of a scaled aircraft model in the Aeronautical and Maritime Research Laboratory 
(AMRL) low-speed wind tunnel. Additional power-off static data were obtained from a 
computational fluid dynamic model for a limited number of aircraft configurations. 

Flight test data for the PC 9/A was required to both validate the flight dynamic model 
and to provide dynamic derivative estimates of the aircraft. A flight test program was 
conducted on a PC 9/A aircraft, serial number A23-045, by the RAAF Aircraft Research 
and Development Unit (ARDU) between November 1998 and February 1999. The details 
of this flight test program and the aircraft instrumentation axe reported in reference [1]. 
The estimation of the stability and control derivatives from the flight test data involved the 
use of system identification techniques, and required control input, aircraft response and 
flight condition data, measured using a high fidelity instrumentation system. Maximum 
likelihood and stepwise regression techniques were employed for the system identification. 
The longitudinal derivatives for the PC 9/A aircraft in the cruise configuration determined 
using system identification techniques are reported in reference [2]. 

This report presents the lateral stability and control derivatives of the PC 9/A extracted 
from flight test data for the cruise configuration, including a comparision with AMRL wind 
tunnel test data and empirical results. Sections 2 and 3 present details of the test aircraft 
and instrumentation system. Section 4 presents the system identification techniques, while 
section 5 discusses the results. 


2 PC 9/A Test Aircraft 

2.1 Aircraft Description 

The flight test aircraft, A23-045 is a Pilatus PC 9/A operated by ARDU. The PC 9/A is 
a single-engine, metal-skinned, low-wing, tandem two-seat training aircraft. The aircraft 
is powered by a Pratt and Whitney PT6A-62 turbo-prop engine flat rated to 950 SHP [3], 
which drives a Hartzell HC-D4N-2A four-blade variable pitch propeller. The aircraft was 
instrumented as outlined in section 3. 
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2.2 Flight Control System 

The aircraft primary flight controls consist of the ailerons, rudder and elevator. The control 
surfaces are manuall y operated from a conventional dual control column and rudder pedal 
arrangement. The stick and rudder pedals are connected to the control surfaces through 
a system of control rods, bellcranks, cables and levers. Trimming control is provided on 
all three axes. 


2.3 Weight, Centre-of-Gravity and Mass Moments-of-Inertia 

The test aircraft, A23-045, was weighed by ARDU prior to the flight tests and had a basic 
mass of 1784.5 kg and a longitudinal centre-of-gravity position of 26.25% MAC [4]. During 
the flight test program, the aircraft weight, centre-of-gravity and mass moments-of-inertia 
varied with fuel usage. These parameters were calculated for each test manoeuvre using 
the equations in Appendix A. 


3 Instrumentation 

Aircraft A23-045 was fitted with an instrumentation system designed specifically for the 
gathering of flight dynamic data. A summary of the instrumentation used in the flight 
test program is included below and a more comprehensive description of the design re¬ 
quirements and calibration is included in references [1] and [5]. The ARDU General Data 
Acquisition System (GDAS) was fitted in place of the rear ejection seat. The data were 
encoded by a 16-bit pulse code modulation system and were recorded onboard the air¬ 
craft using a MARS-2000 14-track tape. Real-time flight test monitoring was provided by 
telemetry data transmitted to the ARDU Primary Analysis Processor (PAP) hut. 

The angular rates (p, q, r) and linear accelerations (ax, ay, &z) were measured using the 
ARDU KAISG1134-1 Motion Platform. This consisted of three Smith Industries 950 RGS 
angular rate gyros and three SunStrand QA1400 servo accelerometers. The aircraft roll, 
pitch and yaw attitude angles (cp, 9 , ip) were obtained by tapping output from the existing 
LISA 2000A Artificial Horizon Reference System (AHRS). 

Static outside air temperature (OAT), indicated airspeed (IAS), angle-of-attack (a) and 
angle-of-sideslip (/3) were obtained from the Rosemount Model 92AN flight test air data 
boom, mounted on the outboard hardpoint on the starboard wing (see figure 1). Aileron 
and elevator deflections were measured using Space Age Control Inc. series 160 cable 
position transducers. Rudder deflection was measured using a type 26V-11CX4C position 
transducer. All sensors were calibrated prior to commencement of the flight test program. 

The aircraft was also instrumented to measure engine torque, propeller speed, inlet turbine 
temperature (ITT), gas generator speed (NG), fuel flow and fuel quantity. The torque, 
propeller speed, altitude and true airspeed, when used in conjunction with a performance 
map for the Hartzell HC-D4N-2A propeller [6], allowed the calculation of engine thrust. 
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4 Methods of Analysis 

During the flight test program, doublet and 3-2-1-1 manoeuvres were performed about 
steady flight conditions at airspeeds between 90 and 200 KIAS. Reference [1] describes 
these manoeuvres in detail. The 26 lateral doublet and 129 lateral 3-2-1-1 manoeuvres were 
analysed using stepwise regression [7] and maximum likelihood estimation techniques [8] 
to determine the lateral stability and control derivatives of the aircraft. 

A range of independant aircraft derivative data were used for comparision purposes and, for 
the maximum liklihood analysis, as a priori estimates to increase the rate of convergence 
of the system identification. This data was souced from AMRL power-off wind tunnel 
tests [9] and empirical data estimated by the University of Sydney [10] using Datcom [11] 
and Pilatus [12] using both Digital Datcom [13] [14] and other ‘long-hand’ methods. 

A right handed orthogonal axes system was adopted for the analysis of the flight test data. 
Positive control surface deflections were defined as elevator trailing edge down, rudder 
trailing edge to port, and starboard aileron trailing edge down, port aileron trailing edge 
up. 

4.1 Stepwise Regression 

Stepwise regression is an unbiased least squares estimator in which new independent vari¬ 
ables are inserted into a model, one at a time, until the regression equation is deemed 
acceptable. The appropriateness of the model can be determined by ex amining a num¬ 
ber of quantities including the squared multiple correlation coefficient and the F statistic. 
The squared multiple correlation coefficient, R 2 , gives a measure of the importance of 
each variable as it is inserted into the equation [7]; however, the improvement in R 2 due 
to the addition of new terms must have some real significance besides simply reflecting the 
inclusion of more terms. This can be determined by monitoring the F statistic, the ratio of 
the regression mean square to the residual mean square. The inclusion of any significant 
terms is generally accompanied by an increase in the F statistic and the best fit with the 
least number of parameters may be obtained by maximising F. Any variable which does 
not make a significant contribution is removed from the model, with the selection process 
continuing until no new variables remain to be inserted into the equation. Whilst stepwise 
regression gives estimates of the derivatives included in the regression equation, it also 
permits a suitable structure for the model to be determined. In addition, it provides an 
independent check on the data estimated using the maximum likelihood technique. 

The stepwise regression technique was applied using code available in the MATLAB© 
Statistics Toolbox [15]. The following model equations were identified during the analysis. 


C Y = Cy 0 + Cy p /3 + Cy St S r 


(1) 


Cn = C„„ + C n ,0 + + C„ r ^ + C„,J r 


( 2 ) 
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Ci = C l0 + Q s /3 ■ 


Cl *W + Clr 2V +ClSa 


6 a 


(3) 


4.1.1 Error Band 

Included in the stepwise regression analysis is the calculation of an error band on the 
estimated derivatives. For a confidence interval of 95%, this error band is approximately 
equal to two standard deviations. Figures 2 to 21 show the stepwise regression derivative 
estimates, including the calculated error band. 


4.2 Maximum Likelihood 

The maximum likelihood estimation technique was applied using the computer program, 
pEst, developed at NASA Dryden Flight Research Center [16]. pEst is an interactive pa¬ 
rameter estimation program which solves a vector set of time-varying, ordinary differential 
equations of motion. 

Prior to the application of the maximum likelihood technique, stepwise regression was 
used to determine a suitable structure for the lateral model equations. Those derivatives 
that were found to be statistically insignificant were not estimated using pEst. However, 
nnlilcp Stepwise regression, these derivatives were not removed from the model equations 
but were instead fixed at their a priori values. C Yp was fixed at a value of -0.16 (per 
radian), C Yr at 0.34 (per radian), C Yga at 0.0 (per degree), C nSa at -0.000035 (per degree) 
and Ci 5 fixed at 0.00035 (per degree). These a priori values were obtained from empirical 
data of references [10] and [12]. The lateral derivatives which were identified using the 
maximum likelihood technique are given in table 1. 



Side Force 

Yawing Moment 

Robing Moment 

Aerodynamic 

Cy 0 

Cn 0 

c h 


- 

C np 

C h 


- 

C nr 

a 

Control 

- 

- 

Ci So 


C Y St 

Cn*. 

- 


Table 1: Lateral derivatives estimated 

The resulting total force and moment coefficient equations used in the state and response 
equations for the maximum likelihood analysis of the PC 9/A axe given in equations 4, 5 
and 6. The differences between equations 1, 2 and 3, and 4, 5 and 6 are important, and 
arise from the inclusion of the a priori values in the maximum likelihood technique. 


pb 


rb 


C Y = C Yo + C Y0 /3 - 0.16^ + 0.34— + C y , r <5, 


2V 


(4) 


4 
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Cn = C no +C ng /3 + C np jy + Cnr^y- 0.000035<f a + C n , r Sr (5) 

Cl = c l0 +C h 0 + C lp ^ + Ci r ~ + C h J a + 0.00035<5 r (6) 


4.2.1 Cramer-Rao Bounds 

For the estimated parameters, pEst calculates a measure of the estimation certainty known 
as the Cramer-Rao bound. A detailed interpretation of this quantity is given in [17]. The 
Cramer-Rao bounds are shown for each derivative estimated by pEst in figures 2 to 21. 
The Cramer-Rao bounds have been factored in accordance with the procedures described 
in [17] to account for the presence of band-limited noise. 


5 Results and Discussion 

Lateral derivatives estimated via maximum likelihood parameter estimation, stepwise re¬ 
gression, empirical methods and power-off wind tunnel experiments are summarised in 
figures 2 to 21. 


5.1 Angle-of-Sideslip Derivatives 

Figures 2 to 6 show the angle-of-sideslip derivatives Cy $ , C n0 and Ci 0 plotted against 
a. The flight test estimates of the side force due to angle-of-sideslip derivative, Cy g , 
generally fall between the empirical and wind tunnel estimates with the results showing 
a small degree of scatter at angles-of-attack above 6°. The estimates of Cy 0 from the 
stepwise regression method are slightly smaller in magnitude than those from the maximum 
likelihood technique. 

Estimates of the directional stability derivative, C n$ , show good agreement between results 
from stepwise regression and maximum likelihood techniquess for both 3-2-1-1 and doublet 
manoeuvres. However, at low angles-of-attack, the flight test estimates are smaller than 
the wind tunnel and empirical estimates, which are both power-off. This difference is 
due to the imbalance between two competing effects acting on the portion of the vertical 
tail within the propeller slipstream, these being the increased dynamic pressure and the 
decreased angle of incidence. The reduction in angle of incidence dominates, resulting in an 
overall reduction in the directional stability of the aircraft. The flight test results show less 
variation with angle-of-attack than the wind tunnel results because the natural reduction 
in directional stability with angle-of-attack is partially balanced by the increasing portion 
of the vertical tail which is subject to increased dynamic pressure within the slipstream. 
In general, the Cramer-Rao bounds and stepwise regression error bands for C n0 estimates 
are small. 
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The maximum likelihood and stepwise regression estimates of the rolling moment due to 
angle-of-sideslip derivative, C t$ , show good agreement with the empirical and wind tunnel 
estimates for both the 3-2-1-1 and doublet manoeuvres. The flight test estimates of C if3 
show little variation with angle-of-attack, as predicted by the wind tunnel results. Similarly 
to Cy s estimates, C ip estimates from the stepwise regression technique are slightly smaller 
in magnitude than those obtained by the ma ximu m likelihood technique. 


5.2 Yaw Rate Derivatives 

Figures 7 to 11 show the yaw rate derivatives, C Ut and Q r , plotted against a. Angular 
rates were not measured during the wind tunnel test program and hence no wind tunnel 
yaw rate derivatives are available for comparision. As mentioned previously, Cy r was 
constrained to its a priori value for the maximum likelihood analysis, and not estimated 
in the stepwise regression analysis. Investigations were undertaken to assess the relative 
contribution of each derivative to the aircraft’s forces and moments to ensure that those 
set to their a priori values were small contributors. This analysis, which is detailed in 
Appendix B, showed the yaw rate contribution to side force to be small. 

The flight test estimates of the yawing moment due to yaw rate, C nr , agree well with the 
empirical estimate for both the 3-2-1-1 and doublet manoeuvres . A high degree of scatter 
is shown in the derivative estimates above 6° angle-of-attack. The stepwise regression 
technique gives estimates of C Ur which are slightly larger in magnitude than the estimates 
from the maximum likelihood technique. 

The rolling moment due to yaw rate, Ci r , whilst being small is an important lateral- 
directional cross-coupling derivative. The flight test anaylsis results show C ir estimates to 
be generally larger than the empirical estimate for both 3-2-1-1 and doublet manoeuvres at 
all angles-of-attack. The flight test estimates of C ir show an increase with angle-of-attack. 
This is an expected result caused by the increase in lift coefficient with increasing angle- 
of-attack, which results in a greater rolling moment as the aircraft yaws (refer [18]). The 
flight test results show a small degree of scatter, particularly for the estimates obtained 
by the maximum likelihood technique. 


5.3 Roll Rate Derivatives 


The roll rate derivative estimates, C Up and Ci p , are plotted against a in figures 12 to 16. 
The side force due to roll rate, Cy p , is a small derivative, not estimated in the stepwise 
regression or maximum likelihood analysis. 

The yawing moment due to roll rate is an important lateral-directional cross-coupling 
derivative. The maximum likelihood estimates of C Up agree well with the stepwise regres¬ 
sion estimates, and estimates from both techniques have small uncertainty levels. However, 
the flight test estimates are larger in magnitude than the empirical estimate. Similarly to 
C lr , C np is influenced by the angle-of-attack and lift coefficient and its value is expected to 
vary with angle-of-attack [18]. This trend is seen in the flight test estimates, which show 
the value of C np to decrease with angle-of-attack. The results also show a small degree of 
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scatter at angles-of-attack above 6°. 

The roll damping derivative, Ci p , estimates from both the stepwise regression and maxi¬ 
mum likelihood techniques show good agreement with the empirical estimate. As expect¬ 
ed [18], the Ci p estimates do not show any variation with angle-of-attack. The results 
display little scatter with small error bands, indicating a high level of confidence in the 
results. 


5.4 Control Derivative 

The control derivatives estimates, Cy 5i .,CW and Ci Sa , are plotted against a in figures 17 
to 21. The yawing moment due to aileron derivative, C nSa , and rolling moment due to 
rudder derivative, Ci Sr , are both small derivatives, which were not estimated in either the 
stepwise regression or maximum likehood analysis. 

Estimates of the side force due to rudder deflection derivative, Cy Sr , are scattered about 
the empirical and wind tunnel estimates. The maximum likelihood estimates of Cy s have 
large Cramer-Rao bounds indicating a higher level of uncertainty in the estimates than 
other derivatives. However, Cy Sr has only a small contribution to the total aircraft side 
force coefficient (refer Appendix B) and hence, the increased uncertainty in this derivative 
will not significantly affect the overall aircraft model characteristics. 

Estimates of the rudder control effectiveness, C nSr , show good agreement with the empir¬ 
ical and wind tunnel estimates and small uncertainty levels, particularly at low angles-of- 
attack. The flight test estimates of the aileron control effectiveness derivative, Ci Sa , fall 
between the empirical and wind tunnel estimates for the 3-2-1-1 and doublet manoeuvres. 
These results exhibit only a small degree of scatter and show a good correlation between 
the stepwise regression and maximum likelihood estimation techniques. 


5.5 Variation of Lateral Derivatives with Rate of Climb and 
Angle-of-Attack 

The variation of lateral derivatives with angle-of-attack and nominal rate of climb for a 
selection of 3-2-1-1 manoeuvres is shown in figure 22. A linear curve fit was applied to 
each rate of climb data set to highlight the trends shown by the data. The influence of rate 
of climb with increasing angle-of-attack is seen in figure 22, with C n&r , Cy g and C n$ all 
increasing in magnitude. These trends are due to increases in the aircraft thrust coefficient 
and the dynamic pressure in the propeller slipstream as the rate of climb increases for a 
given angle-of-attack. These trends are supported by previous theoretical and experimen¬ 
tal studies detailed in references [19] and [20]. For the manoeuvres in figure 22, thrust 
coefficient varies from approximately 0.02 to approximately 0.12. 
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5.6 Manoeuvre Type and Altitude Dependancy 

Tables 2 and 3 show the mean values across the linear angle-of-attack range (0° - 6 ) of the 
aerodynamic derivatives estimated by both maximum likelihood and stepwise regression 
techniques, for each manoeuvre type and altitude band. 


■ 

1 

3-2-1-1 

10 000 ft 

3-2-1-1 

15 000 ft 

Doublet 
5000 ft 

Doublet 

15 000 ft 

wsm 

-1.35E-02 

-1.35E-02 

-1.40E-02 

-1.40E-02 

-1.38E-02 

Cy s 

3.29E-03 

3.00E-03 

3.26E-03 

3.27E-03 

3.63E-03 


1.41E-03 

1.39E-03 

1.33E-03 

1.44E-03 

1.42E-03 

,l P 

C np 

-6.32E-02 

-6.29E-02 

-5.28E-02 

-5.40E-02 

-6.06E-02 

C nr 

-2.01E-01 

-1.97E-01 

-2.03E-01 

-1.80E-01 

-1.91E-01 

Cn s 

-2.02E-03 

-1.91E-03 

-1.89E-03 

-1.73E-03 

-1.86E-03 

C h 

-1.52E-03 

-1.55E-03 

-1.49E-03 

-1.37E-03 

-1.36E-03 

c. 

-5.08E-01 

-5.04E-01 

-4.94E-01 

-4.64E-01 

-4.77E-01 

*p 

Cl, 

1.14E-01 

1.20E-01 

0.962E-01 

1.31E-01 

1.52E-01 


-1.89E-03 

-1.91E-03 

-1.79E-03 

-1.48E-03 

-1.71E-03 


Table 2: Summary of derivatives estimated using maximum likelihood 



1 

3-2-1-1 

10 000 ft 

3-2-1-1 

15 000 ft 

Doublet 
5000 ft 

Doublet 

15 000 ft 

Cy 3 

-1.26E-02 

-1.32E-02 

-1.24E-02 

-1.25E-02 

-1.27E-02 

Cy s 

3.41E-03 

3.50E-03 

3.34E-03 

3.45E-03 

3.68E-03 

Cn B 

1.45E-03 

1.47E-03 

1.45E-03 

1.52E-03 

1.50E-03 

C np 

-5.54E-02 

-5.21E-02 

-5.19E-02 

-4.72E-02 

-5.55E-02 

Cflr 

-2.43E-01 

-2.23E-01 

-2.24E-01 

-2.19E-01 

-2.28E-01 

Cn s 

-1.99E-03 

-1.91E-03 

-1.91E-03 

-1.87E-03 

-1.88E-03 

,t 'or 

C h 

-1.38E-03 

-1.28E-03 

-1.37E-03 

-1.34E-03 

-1.30E-03 

*P 

Ci 

-4.89E-01 

-4.44E-01 

-4.66E-01 

-4.38E-01 

-4.61E-01 

ip 

Ci r 

1.43E-01 

1.43E-01 

1.57E-01 

1.47E-01 

1.63E-01 


-1.82E-03 

-1.65E-03 

-1.73E-03 

-1.61E-03 

-1.70E-03 


Table 3: Summary of derivatives estimated by stepwise regression 

The flight test derivative estimates from the 3-2-1-1 and doublet manoeuvres show similar 
trends for all lateral derivatives. In general, there was no significant change in uncer¬ 
tainty levels, derivative values or degree of scatter in the results for the two manoeuvre 
types at angles-of-attack below 6°. As no doublet manoeuvres were performed above this 
angle-of-attack, no conclusions can be draw as to the influence of manoeuvre type on 
derivative estimates at higher angles-of-attack. As expected, figures 2 to 21 do not show 
any significant altitude dependancy as the results are non-dimensionalised by dynamic 

pressure. 
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6 Conclusions 

The lateral stability and control derivatives of the PC 9/A in the cruise configuration 
have been determined from flight test data measurements using maximum likelihood and 
stepwise regression estimation techniques. Comparisions have been made with derivative 
estimates from empirical methods and wind-tunnel experiments. 

The sideslip derivatives, Cy e and Ci 3 , showed good agreement between the flight test 
and ground-based estimates, while the results for showed small differences due to the 
exclusion of power effects in the wind tunnel and empirical estimates. Flight test estimates 
of the yaw rate derivatives, C nr and Q r , showed good results when compared to empirical 
estimates, and displayed the expected trends with angle-of-attack. Estimates of C Up were 
slightly larger in magnitude than the empirical estimate, while flight test and empirical 
estimates of Ci p compared well. 

The control derivatives, Cy Sr , C nSr and C[ Sa , were also estimated from the PC 9/A flight 
test measurements. Cn Sr and Ci Sa derivative estimates compared well with results esti¬ 
mated from empirical methods and wind tunnel experiments, while Cy Sr estimates showed 
a higher level of uncertainty than other derivatives. 

Trends were found in the Cn Sr , Cy g and Cnp estimates when plotted as a function of 
airspeed and rate of climb that agree with predictions from previous theoretical and ex¬ 
perimental studies. In general, the flight test results did not show any significant variation 
with altitude or manoeuvre type. 
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Figure 1: Air data boom installation. 
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O Flight Test (pEst estimates) 

3 Flight Test (swr estimates) 

.Empirical estimate 

-Wind tunnel estimates (Power off) 





Figure 3: PC 9/A angle-of-sideslip derivatives, 3-2-1-1 manoeuvres at 10 000 ft. 
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Figure 5: PC 9/A angle-of-sidslip derivatives, doublet manoeuvres at 5000 ft. 
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Figure 7: PC 9/A yaw rate derivatives, 3-2-1-1 manoeuvres at 5000 ft. 
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Figure 9: PC 9/A yaw rate derivatives, 3-2-1-1 manoeuvres at 15 000 ft. 
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O Flight Test (pEst estimates) 
□ Flight Test (swr estimates) 

- Empirical estimate 




Figure 11: PC 9/A yaw rate derivatives, doublet manoeuvres at 15 000 ft. 
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Figure 15: PC 9/A roll rate derivatives, doublet manoeuvres at 5000 ft. 
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Figure 16: PC 9/A roll rate derivatives, doublet manoeuvres at 15 000 ft. 
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Figure 17: PC 9/A control derivatives, 3-2-1-1 manoeuvres at 5000 ft. 
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Figure 19: PC 9/A control derivatives, 3-2-1-1 manoeuvres at 15 000 ft. 
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Figure 21: PC 9/A control derivatives, doublet manoeuvres at 15 000 ft. 
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Appendix A: Weight, Centre-of-Gravity and 
Mass Moments-of-Inertia 


The aircraft mass, centre-of-gravity and mass moments-of-inertia were determined for each 
test man oeuvre based on the aircraft fuel mass. The basic aircraft, pilot and fuel masses 
and moment arms from references [3] and [4] are given in table Al. The moment arms are 
measured relative to the aircraft datum, located 3m forward of the engine firewall, and 
2m below the forward fuselage reference plane. 



Mass (kg) 

Moment Arm (m) 

Basic Aircraft 

1784.5 

4.306 

Pilot 

81.6 

4.061 

Fuel 

Variable 

4.178 


Table Al: Flight test aircraft mass distribution [3, 4]- 

The mass of the aircraft is calculated as the sum of the basic aircraft, pilot and remaining 
fuel masses, by the following equation. 


M = 1866.1 + M f (kg) (Al) 

The longitudinal centre-of-gravity position of the flight test aircraft was determined from 
the foDowing equation. 


x c.g. 


_ 8015.5+4.178Mf / \ 
~ 1866.1+M F ^ 


(A2) 


During the flight test program, the aircraft longitudinal centre-of-gravity position varied 
between 4.27 m (24.2%MAC) and 4.30 m (25.6%MAC). The lateral and vertical centre- 
of-gravity positions were assumed to be invariant with fuel usage, and were fixed at values 
of 0.024 m and -2.2 m, respectively, relative to the aircraft datum. 

The mass moments-of-inertia of the aircraft are given by the following equations. 


Ixx 

= 2505.9 -1- 6.177 (M - 

- 1866.1) 

(kg.m 2 ) 

(A3) 

Iyy 

= 6622.2 + 0.033 (M - 

- 1866.1) 

(kg.m 2 ) 

(A4) 

Izz 

= 8467.1 + 6.188 (M - 

-1866.1) 

(kg.m 2 ) 

(A5) 
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Ixy = 49.0 + 0.0057 (M - 1866.1) (kg.m 2 ) 


I xz = 196.9 + 0.0041 (M - 1866.1) (kg.m 2 ) 


I YZ = 3.0 + 0.0007 (M- 1866.1) (kg.m 2 ) 

> ' 


(A6) 


(A7) 


(A8) 
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Appendix B: Side Force, Rolling Moment and 
Yawing Moment Coefficients 

An analysis was carried out on the PC 9/A flight test data to determine the relative 
contribution of each lateral derivative to the aircraft’s force and moment coefficients, C Y , 
C n and Q. This analysis was undertaken to confirm that contributions from derivatives 
set to their a priori values in the maximum likelihood analysis were not significant in the 
total force and moment equations of the aircraft. This analysis also confirmed the stepwise 
regression results, which showed several derivatives to be insignificant to the regression 
equation. 


B.l Side Force Coefficient 

The side force coefficient may be represented by: 


C Y = Cy 0 + C Yg P 4- C Yp + C Yr + C Ysa 6 a + C Ygr 6 r (Bl) 

A comparision of the relative size of each component of the above equation is shown in 
figure Bl for a 3-2-1-1 manoeuvre. The results show that changes in side force are mainly 
due to the change in the aircraft’s angle-of-sideslip. The roll rate and rudder deflection 
contributions to side force are small, while yaw rate and aileron deflection contributions 
are negligible. 


B.2 Yawing Moment Coefficient 

The yawing moment coefficient may be represented by: 

Cn = Cn o + C ng P + Cn p ^ r— ^ + C Ur (^y^J + ^ n S a + ^ n Sr 

A comparision of the relative size of each component of the above equation is shown in 
figure B2 for a 3-2-1-1 manoeuvre. The results show that the yaw stiffness derivative, C ng , 
and the yaw control derivative, C nSr , axe important, while the remaining derivatives are 
less significant to the total yawing moment coefficient of the aircraft. 
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B.3 Rolling Moment Coefficent 

The rolling moment coefficient may be represented by: 
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C = C h + <V + Q, + c lr + <v. + c w #r (B3) 

ansion of the :;elative size of each component of the above equation is shown in 
-2-1-1 manoeuvre. With the exception of Ci 5r , each of the rolling moment 
'cant to the total rolling moment coefficient of the aircraft. 



Figure Bl: Relative magnitudes of Cy equation components. 
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Figure B2: Relative magnitudes of C n equation components. 



Figure B3: Relative magnitudes of Ci equation components. 








DISTRIBUTION LIST 


A Correlation between Flight-determined Lateral Derivatives and Ground-based Data for 
the Pilatus PC 9/A Training Aircraft in Cruise Configuration 

Hilary A. Keating, Nick van Bronswijk, Andrew D. Snowden and Jan S. Drobik 


DEFENCE ORGANISATION 


Number of Copies 


Task Sponsor 
COPS HQAC 
S&T Program 

Chief Defence Scientist 'j 

FAS Science Policy > 

AS Science Corporate Management J 
Director General Science Policy Development 
Counsellor, Defence Science, London 
Counsellor, Defence Science, Washington 
Scientific Adviser to MRDC Thailand 
Scientific Adviser Policy and Command 
Navy Scientific Adviser 
Scientific Adviser, Army 
Air Force Scientific Adviser 
Director Trials 

Aeronautical and Maritime Research Laboratory 

Director, Aeronautical and Maritime Research Laboratory 
Air Operations Division 

Chief, Air Operations Division 

Research Leader, Avionics and Flight Mechanics 

Head, Flight Mechanics Applications 

Task Manager - G. J. Brian 

Author - H. A. Keating 

Author - N. van Bronswijk 

Author - A. D. Snowden 

Author - J. S. Drobik 

B. A. Woodyatt 

K. L. Bramley 

DSTO Research Library and Archives 
Library Fishermans Bend 


1 

1 

1 

Doc Data Sht 
Doc Data Sht 
Doc Data Sht 
1 

Doc Data Sht 
Doc Data Sht 
1 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



Library Maribyrnong 1 

Library Salisbury 2 

Library, MOD, Pyrmont Doc Data Sht 

US Defense Technical Information Center 2 

UK Defence Research Information Centre 2 

Canada Defence Scientific Information Service 1 

NZ Defence Information Centre 1 

National Library of Australia 1 

Capability Systems Staff 

Director General Maritime Development Doc Data Sht 

Director General C3I Development Doc Data Sht 

Director General Aerospace Development 1 

Aircraft Research and Development Unit 

MSSUP1 1 

Intelligence Program 

DGSTA Defence Intelligence Organisation 1 

M anag er, Information Centre, Defence Intelligence Organisa- 1 

tion 

Corporate Support Program 

Officer in Charge, TRS, Defence Regional Library, Canberra 1 

UNIVERSITIES AND COLLEGES 

Australian Defence Force Academy Library 1 

Head of Aerospace and Mechanical Engineering, ADFA 1 

Deakin University Library, Serials Section (M List) 1 

Librarian, Flinders University 1 

Head of Aeronautical Engineering, University of Sydney 1 

Head of Centre of Expertise in Aerodynamic Loading, RMIT 1 

Adrian Pearce, Computer Science, Curtin University WA 1 

OTHER ORGANISATIONS 

NASA (Canberra) 1 

Info Australia 1 

ABSTRACTING AND INFORMATION ORGANISATIONS 

Library, Chemical Abstracts Reference Service 1 

Engineering Societies Library, US 1 




Materials Information, Cambridge Science Abstracts, US 1 

Documents Librarian, The Center for Research Libraries, US 1 

INFORMATION EXCHANGE AGREEMENT PARTNERS 

Acquisitions Unit, Science Reference and Information Service, 1 

UK 

Library - Exchange Desk, National Institute of Standards and 1 

Technology, US 

National Aerospace Library, Japan 1 

National Aerospace Library, Netherlands 1 

SPARES 

H. A. Keating 5 

Total number of copies: 55 



Page classification: UN CLASSIFIED 


DEFENCE SCIENCE AND TECHNOLOGY ORGANISATION 
DOCUMENT CONTROL DATA 


1. CAVEAT/PRIVACY MARKING 


2. TITLE 3. SECURIT 

A Correlation between Flight-determined Later- Document 
al Derivatives and Ground-based Data for the Title 
Pilatus PC 9/A Tr aining Aircraft in Cruise Abstract 
Configuration _ 


4. AUTHORS 

Hilary A. Keating, Nick van Bronswijk, Andrew 
D. Snowden and Jan S. Drobik 


3. SECURITY CLASSIFICATION 


6a. DSTO NUMBER 

DSTO-TR-0988 


8. FILE NUMBER 

Ml/9/746 


5. CORPORATE AUTHOR 

Aeronautical and Maritime Research Laboratory 
506 Lorimer St, 

Fishermans Bend, Victoria, Australia 3207 


6c. TYPE OF REPORT 
Technical Report 


10. SPONSOR 

COPS HQAC 


14. RELEASE AUTHORITY 


11. No OF PAGES 

38 


7. DOCUMENT DATE 
June, 2000 


12. No OF REFS 
20 


Chief, Air Operations Division 


6b. AR NUMBER 

AR-011-475 


9. TASK NUMBER 

AIR 97/135 


13. URL OF ELECTRONIC VERSION 

http: //www.dsto.defence.gov.au/corporate/ 

reports/DSTO-TR-0988.pdf_ 


15. SECONDARY RELEASE STATEMENT OF THIS DOCUMENT 
Approved, For Public Release 

OVERSEAS ENQUIRIES OUTSIDE STATED LIMITATIONS SHOULD BE REFERRED THROUGH DOCUMENT EXCHANGE, PO BOX 1500, 
SALISBURY, SOUTH AUSTRALIA 5108 _ 


16. DELIBERATE ANNOUNCEMENT 

No Limitations 


17. CITATION IN OTHER DOCUMENTS 

No Limitations 

18. DEFTEST DESCRIPTORS 


Aerodynamics 

Aerodynamic configurations 

Aerodynamic stability 

Flight tests 

Lateral stability 

PC 9/A aircraft 

System identification 



19. ABSTRACT 

A series of flight tests were conducted on the PC 9/A aircraft, A23-045, at the Royal Austr alian Air 
Force’s Aircraft Research and Development Unit. System identification techniques were applied to the 
data obtained from these flight tests to determine the stability and control derivatives of the aircraft. 
The lateral results for the aircraft in cruise configuration are presented in this report and comparisons 
are made with empirical and ground based estimates. _ 






















